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Abstract 

A non-perturbative s-wave renormalization of the pion in a hot and baryon 
rich medium is presented. This approach proceeds via a mapping of the canon- 
ical pion into the axial Noether's charge. The mapping was made dynamical 
in the Hartree-Fock-Bogoliubov random phase approximation (HFB-RPA) . It 
is shown that this approach, while order mixing, is still symmetry conserving 
both in the baryon free and baryon rich sectors, at zero as well as finite tem- 
perature. The systematic character of this approach is emphasized and it is 
particularly argued that it may constitute an interesting alternative for the 
non-perturbative assessment of the nuclear matter saturation properties. 
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I. INTRODUCTION 



The understanding of the nuclear equation of state for hot and dense matter is one of 
the fiercely followed objective in the nuclear problem. Early studies although successful 
in quantitatively assessing the nuclear matter ground-state properties, rely heavily on phe- 
nomenological input, as the nuclear two body interaction for instance [0,0- These bear very 
little to the nowadays admitted fundamental theory of the strong interaction, the quantum 
chromodynamics (QCD). Due to its strong coupling limit in the low energy regimes, QCD 
is in fact intractable as a theory of its fundamental degrees of freedom. On the other hand, 
QCD has also symmetries which are as much fundamental. Their realization through effec- 
tive hadronic theories seems to be a wise compromise and a strategy adopted in the so-called 
quantum hadrodynamics (QHD) (see ^ for a review). 

Of particular interest is chiral symmetry. Its nonlinear realization, through chiral per- 
turbation theory (ChPT) for instance, has proven successful all over the past decade in 
understanding the low energy pion physics. This, however, falls short of answering some 
important questions as how to deal with highly collective structures like bound states or 
resonances and how to fulfill the fundamental unitarity condition. Thus one realizes that 
besides QCD's symmetries, the QHD inherits as well the difficult aspects of the strongly 
interacting fields. Any reasonable handling of these theories requires therefore an adequate 
treatment of the complicated vacuum structure. This is addressed in general through non- 
perturbative methods. The symmetry, however, imposes stringent constraints on the way 
the evaluation of the dynamics is conducted. Therefore one is called to be very selective in 
choosing non-perturbative approaches that help in gathering the dynamics without destroy- 
ing the symmetry. 

This task, although trivial in the case of the coupling constant perturbation (CCP), is 
full of subtleties and requires in fact much more effort from within the strong interaction 
physics community. In this regard, as an educated non-perturbative approach, the concept of 
symmetry conserving dynamical mapping (SCDM) was introduced in The idea is to map 
the original Fock-space, created at the quantization and which supports the CCP approach, 
to some ideal Fock-space selected via a given symmetry conserving mapping (SCM). The 
next step is to make a projection onto the physical Fock-space which is a subspace of the 
ideal Fock-space. The so-called projection is realized by making the SCM dynamicalQ. 

There are various advantages in performing non-perturbative calculations following this 
scheme. The most interesting aspect is certainly the possibility of tracking all kind of trans- 
formations taking place in the Fock-space while gathering the dynamics In particular 
the asymptotic particles, which are highly dressed by the dynamics, have also well defined 
normalizable states in a well defined Fock-space. This is a clear departure from some tra- 
ditional non-perturbative approaches in which one loses this precious insight. Probably the 
difficult task in this picture is to actually find the right symmetry conserving mapping. 

Beside the expansion approach which was sorted-out using Holstein-Primakoff map- 
ping, an other SCM, based on a field-to-current mapping, was presented in That this 



Further precautions need, of course, to be observed in projecting out unphysical states. The 
necessary technics to address this question are well available in the literature 
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second mapping is systematic in handling complicated situations was proven in |^ within 
the chirally invariant SU{3) x SU{3) Gellmann-Levy lagrangian. It will be further con- 
solidated here in considering a baryon rich environment at finite temperature within the 
SU (2) X SU{2) version of the model. It will be also made clear later on that the field-to- 
current mapping is a robust SCDM which is able to handle any complicated dynamics as 
long as the symmetry is linearly realized in the lagrangian. In what follows the field-current 
mapping will be made dynamical in the Hartree-Fock-Bogoliubov-RPA scheme. This was 
the approximation adopted previously in [^,0 . Higher approximations are also possible and 
these will soon be reported upon 

This paper is organized as follows: We briefiy review in section III the effective chiral 
model used all over the paper and spell out the basic rule of the game, namely the preser- 
vation of the lowest chiral ward identity throughout our non-perturbative treatment. This 
will be realized via an HFB-RPA scheme. Therefore, we first introduce in subsection-A a 
Bogoliubov rotation and perform a self-consistent mean field calculation which prepares the 
HFB basis. The RPA fiuctuations build on the HFB basis will be evaluated in subsection-B. 
It will be shown that the Goldstone theorem is fulfilled exactly. In section IV, the finite 
temperature and finite chemical potential extensions of this approach are made. Finally, the 
conclusions are drawn in section V. 



II. THE MODEL 



As a starting point, we recall the lagrangian density of the SU (2) x SU (2) linear cr-model 
with nucleonic degrees of freedom. To realize an explicit chiral symmetry breaking in the 
PCAC sense, one adds the usual linear term in the sigma field. 



1 r 1 
C^ = ij[i jd + g{a + zf7r75)] ^ + - [{d^nf + {d^a^ 
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TT + (T 



4 



TT + cr 



ca 
(1) 



Here ip, a and vr represent the bar nucleon, sigma, and pion fields, respectively. To fix our 
notations for latter use, these are given below, in terms of their respective creation and 
annihilation operators, by the usual plane wave expansions 



dp 



(27r)32Ej 



:5:[«(p,r)e-^P%. + t;(p,r)e^P^rf+ 



■pi) , H^) = J 



dp 



v/(27r)32cup 



e^px^p + e-px5+ (2) 



The on-shell energies are taken to be Up = \/ jji^ + p^ for the bosons, and Ep = p for the 
massless fermions. The definition of the bi-spinors u(p, r) and f (p, r) considered here is the 
one which takes care of this limiting case |^. 

As stated in the introduction, the aim here is to conduct a non-perturbative treatment of 
the model that goes beyond the usual coupling constant perturbation (CCP) as well as the 
semi-classical approaches. These two are, presently, the most extensively used approaches 
in the study of the dynamics of these QCD inspired hadronic models. Before embarking 
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in this, we recall the commonly used classical mean-field approximation. In this case, the 
masses of the different particles are given by 

ml = i^ + \^a)\ ml = i? + ?,\\a)\ ml = -g{a), c = i? {a) + X"" {af (3) 

where {a) denotes the vacuum expectation value of the sigma field. At this level, the 
Goldstone nature of the pion, in the chiral limit (c = 0), is manifest. Indeed, one can also 
verify that the Ward identity 

-B.-(0) = ^ (4) 

is fulfilled. Of course, all higher Ward identities (WI) are also preserved by the CCP approach 
at each order of the perturbation, and thus, at this lowest order too. The present WI holds 
as well in the case of the semi-classical non-perturbative approach a la l/A^-expansionQ. In 
fact, in the large A^-limit, the condensate and the pion mass are obtained as solutions of two 
coupled BCS equations in the Hartree-Bogoliubov approximation (see for instance 0,0). 
The fiuctuating part of the sigma field is found to be not involved at all in building the 
variational ground state for the whole 1/A^-expansion approach. As such the sigma mass 
has an exclusively perturbative character. Here too, the whole hierarchy of WI's, as in 
the case of the CCP, is preserved. In the present paper, however, our ambition is rather 
modest. Our aim is simply to set up a non-perturbative approach, which transcends the two 
approximations cited above, but still preserves in a systematic way the lowest WI in Eq.(|), 
at all conditions of temperature and baryon density. 

A. HARTREE-FOCK-BOGOLIUBOV MEAN-FIELD 

The starting point for our mean field considerations is the following Bogoliubov trans- 
formation for both nucleon and meson operators 

ttp = Un{p)ap - v^{p)a^p ; /?+ = u„{p)h^ - v„{p)h^p - Wq6{^) 

Cpr = UF{p)Cp.,-rVF{p)d^pr ] Dpr = UF{p)dl^ + rVF{p)c_pr (5) 

where u.j,{p), Ua{p), Vn{p), Va{p), upip) and vpip) are even functions of their arguments. To 
account for a finite (a) condensate in the Goldstone phase, the Bogoliubov transformation is 
made inhomogeneous for the sigma field by introducing a c-number shift wq. The canonicity 
of the above transformations is enforced by means of the following conditions : 

U^ripf - V^{pf = 1 , U^ipf - V„{pf = 1 , UpipY + Vpipf = 1 (6) 

Accordingly, the vacuum |0) of the theory at the classical mean field, defined by: a,r|0) = 
b„\0) = c|0) = d\0) = 0, is unitarily transformed to a new one 1$), such that 

= = ^1$) = D\^) = , 



is the number of pion charges 
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and explicitly given in terms of the following squeezed vacuum 



1$) = exp 



|0). 



(7) 

where z.„{p) = arcch[M^(p)], z„{p) = arcch[uo-(p)] and zf{p) = arccos[uir(p)]. The Hamilto- 
nian of the model derived from Eq.(|Tp and written in the quasi-particles basis reads: 



n = no + v[/3o + (3^] + / dpSpip) ^ [C+ Cp, + D;^Dpr] + J dp£^{p) 
+ j dp£^{p)p;Pp + I dpcp{p)J2r [C+D+p, + D_p,Cp, 



+ 



fix : 



+ / rfx : 



(8) 



where a represents the shifted sigma-field: a = a — {a), and (a) is the sigma condensate 
which is given in term of the parameters of the Bogoliubov transformation in Eq.(|^) by 

^ (^0 ) + jbo) ^ K(0)+t;.(0))(wo + w*o) 

v/(2^)'2/i V'(27r)32/i 

The semicolons " :" in the expression of the hamiltonian denote operators normal order- 
ing. Finally, the coefficients Hq, t], £f,tt,ct and CiT'^jr.o- read explicitly 

Ep {uF{pf - VF{pf) - 2g{a)uF{p)vF{p) 

u;p{3v^{pf + v„{pf + 2) 
3A2(<t)2 



Ho = -7 
+ 

+ 



dp 



(27r)3 
7^ + 5/^ + 2hh 



[la + h] + ^Ht^ + - C(O^) 



(9) 



where 7 = 4 is the spin-isospin degeneracy of the nucleon, and 
M,(0)+t;,(0) 



2/i 



^y{a)h + 3A'(a)4 + A'((t)^ + ^l'{a) -c-qIf 



Cf 



(p) = 9{^) [uF{pf - VF{pf) + 2EpUf{p)vf{p) 



C^{p) = UJpU^{p)VT,{p) + 



Ca{p) = UpU„{p)v„{p) 



Mp)+Vn{p)) 

2 2ujp 

3A^ Mp) + va{p)y 
2 2tu„ 



5/^ + + {ay 



h + h + {af 



Sf{p) = Ep {uF{pf - VF{pf) - 2g{a)uF{p)vF{p) 

Mp) + v„{p)y 



8^{p) = ujp (u^ipf + v^ipy'j + A 



2u„ 



8„{p) = Up iu^ip) + v^{p) ) + 3X 



Mp) + Vajpy^ 

2oj„ 



5/^ + + {ay 



h + L + {ay 



(10) 



/tt, la and If are expectation values on the squeezed vacuum of bilinear forms of the pion, 
the sigma and the baryon fields, respectively. They read 



dx ($|7ri(x)7ri(x)|<l>) 

dx (<l>|a(x)a(x)|$) = 
rfx ($|V^(x)V^(x)|$) = 



dp ju^jp) + Vn{p)y 
(27r)3 2LUp 

dp {Ua{p) + Va{p)f 

(27r)3 
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dp 



■2uf{p)vf{p) . 



(11) 



To fully fix the amplitudes Ut,, f^, u^, v^, uf, vf, as well as the value of the sigma con- 
densate (a), we make use of the Rayleigh-Ritz variational principle. Minimizing the ground 
state energy Tig = /{^\^) while maintaining the canonicity condition of the Bogoli- 

ubov transformation Eq.(^, one gets, after a straightforward algebra, the following set of 
equations 



Su^{p) 



6ua{p) 



5uf{p) 



6{a) 










cf{p) 









. 



(12) 



From the above, one sees that at the minimum of the HFB ground state, the Hamiltonian 
in Eq. (|) reduces to a sum of un-coupled quasi-particle modes (bilinear parts) and residual 
interactions which are normal ordered with respect to the squeezed vacuum. The set of 
equations in Eq.(|T^) constitutes four coupled and self-consistent gap equations. Inserting 
the solutions of these into the definitions of the quasi-particle energies, one gets the following: 

S^{p) = Up {u^{p) - v^{p)f = + [5/^ + h + (o")^] , 

8a{p) = UJp (uaip) - Va{p)f = ^cu2 + 3A2 + , 



£f{p) 



Ep 
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Uf{p) -VpipY 



p^ + g'^{cry 



lF + ^\\h + h) + \\(y)^ . 



(13) 



{a) - {ay 

Finally, one deduces the BCS gap equations for the quasi-particle masses as well as the 
condensate 

s^{o) = fi' + x'[5i^ + ia + {ay] , 



^f(O) 
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■If + 



I. 



3A^ X'jay 

— //^ — /i^ 



(14) 
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Here /tt, la and Ip are nothing but the one loop tadpoles of the fully dressed quasi-pion, 
quasi-sigma and quasi-fermion propagators, respectively. 

dp 1 



la 



(27r)3 2^p2 + ^2(0) 

dp 1 
(2vr)3 2 /p2 + ^2(o)^ 



/ dp SpjO) 



The Feynmann diagrammatic representation of the gap equations in Eq.(p!4D is displayed in 
Figure 1. Through the Yukawa coupling, the fermion is in fact minimally coupled to the 
rest of the gap equations. The gap equation for the fermion could be more tightly coupled 
if one considers, for instance, a chirally invariant quartic interaction among fermions a la 
Nambu Jona-Lasinio. This will be briefly discussed at the end of the paper. 

At this point, a comment is in order. After inspection of the BCS gap equations Eq.(|l^), 
it is clear that a solution with finite condensate {{a) ^ 0) will lead to the following values 
for the masses of the quasi-pion and quasi-sigma : 

^^(0) = ^^ + 2A^ [/.-/.] 

^^(0) = ^^ + 2A^(a)^ (16) 

In this Nambu-Goldstone phase, these two modes are not degenerate as it should be. There- 
fore, the difference — la as well as the fermionic tadpole Ip are in fact finite. Thus in 
the chiral limit (c — 0), the Goldstone mode seems to be absent from the spectrum of 
the theory, since Et, does not vanish. Consequently, one can check that the Ward identity 
in Eq.® is also violated. Since the condensate and the Goldstone mass are the building 
blocks for all higher Ward identities, it is clear that the quasi-sigma mass is also unphysical 
in this approach. It appears then that the squeezed vacuum, apart from the fact that it can 
accommodate a finite condensate, does neither have the right curvature for the sigma-like 
excitation, nor a valley along the pion-like excitation. In the following, we can only propose 
a way to resolve the second problem, namely correcting for the pion-like excitation so as 
to fulfill the Goldstone theorem. This will have quantitative consequences, as there will be 
important corrections to the ground state energy. 



B. RPA FLUCTUATIONS 



To proceed further, we recall the Hamiltonian in Eq.@. In the HFB basis, the latter 
takes the form: 



n = no + l dpSpip) J2 [CprCpr + D^^D^,] + J dpS^{p)a^d^ + j dpSa{p)l3p, 



(ix : 



(17) 
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Here the fields 7rj(x) and cr(x) have now the following plane wave expansions in 

the quasi-particle basis 



^ ■ /(27r)32^^(p) r ^ 



dp 

7r,(x) = / (e^P'^ap, + e-^P^«+ ) , 



^ / lira 



where the harmonic modes are the respective quasi-particle energies. 

For a further gathering of the dynamics, we want to proceed by a perturbative diago- 
nalization of the Hamiltonian^. This can now be done via two kinds of RPA's; one with 
sigma-like excitations and an other with pion-like excitations. For symmetry considerations, 
we rather favor the second. Indeed it was shown in that a RPA with the quantum num- 
bers of the pion field can be made symmetry conserving by enlarging the RPA excitation 
operator in such a way to accommodate the whole set of excitations formally present in 
the two-body axial charge operator. On the other hand, Noether theorem defines the axial 
current to be 

At, = 1^1,1 J-^i^ + d.aix- - d,Ti^{a + (a)) . (19) 

The symmetry operator Q5, being the volume integral of the time component Ajq, takes 
in second quantization a structure in which one can identify linear forms standing for the 
creation and annihilation of quasi-pion modes at rest, as well as bilinear forms corresponding 
to creation and annihilation of pairs of either fermions or bosons with opposite parities and 
vanishing total momentum. 

According to the Goldstone theorem, the action of the axial charge on the full correlated 
vacuum with a spontaneously broken phase creates in the chiral limit a non normalizable 
state which is nothing but the pion state 

Ql\vac) oc Ivr'^) . 

Here one can mimic this situation and build, on an approximate vacuum which is considered 
here to be a RPA ground state {\RPA)), a normalizable pion state, away from the chiral 
limit, using the following RPA excitation operator 

k k 

+ J2[xl'\k)ClDt,r'-Yj'\k)D_,r'C,r] + j:[xl'\k)ClC,r' -Yj'\k)D,r'Dt] . (20) 

krr' krr' 



^What is meant here is not the coupling constant perturbation. This word is used here to highUght 
the additive character of this approach as opposed to a variational self-consistent approach. 
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The latter is defined to contain the same pair-excitations as the axial charge but with 
amplitudes which remains to be fixed dynamically. The RPA ground state is defined as 
the vacuum for the above operator Q'^\RPA) = 0, while the single pion state, given by 
Q'^~^\RPA) = |7r°), is normalized as follows 



{RPA\ fQ" \RPA) = 6, 



(21) 



Using Rowe's equation of motion method |Tl]|, one is lead to the usual RPA secular 
equations 



{RPA\ \6Ql, \H,Qt+]] \RPA) = u,{RPA\ \6Ql,Ql+] \RPA) . 



(22) 



These are solved within the usual quasi-boson approximation in which the \RPA) ground 
state is replaced by the HFB ground state |$). This procedure which seems brutal and 
inconsistent is nevertheless commonly used and fully under control. It can be shown that it 
is equivalent to a lowest order truncated bosonization of all the bilinear operator products 
which form the Hamiltonian and the excitation operator. In the ansatz above, only the 
number non-conserving bilinear terms contribute effectively to the RPA equation Eq.(^). 
The remaining number conserving terms decouples totally. They should appear in higher 
non harmonic orders in the bosonization. This point will not be developed further here. The 
equations of motion Eq.(p2D are in fact a set of coupled channel equations where all possible 
states with opposite parities scatter according to the vertices allowed by the Lagrangian 
Eq.(|l]). These coupled-channel Lippmann-Schwinger equations are further coupled to a 
Dyson equation for the pion propagator. These processes are represented schematically 
below 



TT ^ TT 
TTO" — S> TT 



TTa — > TTa 




TT 







Since the linear cr-model does not provide for a coupling between pairs of fermions and 
bosons the matrix which governs the subspace for pure fermionic re-scattering is actually 
diagonal. Including quartic interaction terms a la Nambu Jona-Lasinio leads to a fermionic 
RPA type of re-scattering. This is briefly shown at the end of the paper. 

Let us turn now to the proper resolution of the RPA equations. These can be recast in 
the following eigenvalue problem form 



A{qr,q,) \VM) 




(23) 



As anticipated above, the amplitudes X'^^\ Y^^\ X^^^ and Y^^^ in Eq.(^), corresponding to 
the excitations generated by all number conserving pairs of operators, are found to decouple 
in the present case. It will be seen later on that this situation does not persist at flnite 
temperature for the number conserving bosonic pairs. The number conserving fermionic 
pairs, on the other hand, start contributing only at flnite total three momentum and flnite 
baryon density. Thus in the eigenvalue RPA equation above, A and B are 3x3 matrices. 
The A matrix is the sum of a diagonal matrix which represents the free propagation of 
re-scattering states, and the B matrix such that 
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S^{0) 

AquQ2)=\ [£M)+£M)] 

2£Fiqi] 



(24) 



The interaction is fully encoded in the B matrix which reads 



B{qi,q2) 



2^.(0) 



/ 

2A2(a)7^(gl 
V ^(27r)-i 



2A2(a)7^(g2) 
2A2 7^(gl)7^(g2) 




7^(2vr) 





(25) 



with 



-1/2 



(26) 



Finally the norm matrix M as well as the RPA eigenvectors Uu{q) and Vy{q) are given by: 



h 




I V \ 



uAq) = xi^\q) , = K^'Kq) 



\xi'\q)) 



\Y!;'\q)) 



(27) 



where Id is the 3x3 identity matrix. 

The solution of the eigenvalue problem proceeds through straightforward calculation. 
This can be carried analytically here via a so called Feshbach projection from the two-particle 
states subspaces onto the single pion-state subspace. We get the following expression for the 
RPA eigenvalues 



2A2 



(28) 



where T^pp and S^o- are the RPA bubbles given by 

dp A£f{p) 



'FF 



.UJ„ 



[uj: 
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(27r)3 _ 4^^(p)2 

dp £^{p)+£^{p) 



(29) 



(27r)3 2£Ap)£Ap) ^1 - {£^{p) + £„{p)f 

Eqs. ([28|j29| ) are diagrammatically represented in Fig. (2). As it was motivated earlier, one of 
the RPA solutions has to have the Goldstone character. This solution is commonly known, 
in the nuclear problem, as the spurious solution. This negative connotation expresses the 
fact that this solution does not correspond to a real excitation of the system, but rather is an 
expression of a dynamical breaking of a symmetry in the system as, for instance, the space 
rotation or translation of a nucleus. In the present case, we hold the spurious solution for a 
real excitation of the vacuum and it corresponds to the pionic mode. To retrieve explicitly 
this solution, we first notice the following identities 



S..(0) 



£m-^m 



Sff(O) 



£f{Q) 



(30) 
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These, together with the gap equations Eq.(|T4D, allow to rewrite the RPA frequencies in 
Eq.(p8l) into the following form 



+ 



s..(0)] 



1 



2A2 



[UJt. 



(31) 



It is clear, from the expression above, that the Goldstone solution is manifestly present in 
the RPA spectrum. Besides, there exist evidently a continuum of solutions which correspond 
to two cuts; one is the free propagation of pairs of fermions with opposite parities, and the 
other corresponds to the quasi-sigma and quasi-pion scattering process. 

It is well known from the nuclear many-body problem that the random phase approx- 
imation is the perturbative procedure which further diagonalizes the residual interaction 
inherited from the mean field calculation. In the present case, the considered RPA did not 
affect the whole residual interaction. Actually, the part of the Hamiltonian which was diag- 
onalized here is the one which allows the transition between states of mixed parities. This 
part is given by 



+ j dyi: 



)2 



(32) 



and the full Hamiltonian reads 



7i = 7^0 

'Hresidual ' dx 



T^RPA + 'Hresidual 

A2 



A^(a)a^ + -((7fr + a^j-^7^a^ 



(33) 



where Tio is the mean field contribution to the ground state energy and the semicolons stand 
for the normal ordering with respect to the HFB squeezed state. The RPA diagonalization 
brings as well a finite contribution to the vacuum energy. In the RPA basis, the TiRPA piece 
of the Hamiltonian has, therefore, the following form 



'Hrpa 

'Hrpa '■ 



Erpa + : 'Hrpa '■ 



q,u 



Erpa = {RPAIHrpaIRPA) 



(34) 



where the semicolons indicate this time a normal ordering with respect to the RPA ground 
state. The evaluation of the RPA contribution to the ground state requires in fact a finite 
three momentum RPA calculation which was not addressed in this paper. Therefore the 
Erpa will not be explicitly given here. This point as well as the nuclear matter equation 
of state will be addressed soon in a forthcoming publication. As a prerequisite, let us next 
see how to extend the present approach to the case of finite temperature and finite baryon 
density. 
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III. FINITE TEMPERATURE AND BARYON DENSITY EXTENSIONS 



The extension of the above approach to a baryon rich system at finite temperature 
follows through very standard technics. One of these is the Thermo Field Dynamics (TFD) 
[T^ . This approach embraces very well the concept of dynamical mappings since it keeps 



transparent the notion of the vacuum and thus of the Fock space. However, it remains very 
much involved in the present type of calculations because it proceeds via a doubling of the 
dynamical variables {i.e. doubling of the Fock space). In what follows we adopt the more 
practical and rather standard approach which consists in evaluating the Grand Canonical 
Potential Q of the grand canonical ensemble. The conditions of stationarity (vanishing 
of the first derivative of Q), and stability (positivity of its second derivative) lead to the 
thermodynamical equilibrium of the system. In the present case, Q takes the form 

n= {H) - TS - /i(iV) . (35) 

Given the grand canonical partition function Z = Tr[e~'^^^"^'^^], and density operator 
D = Z~^[e~^^^~^^'^] , with P = l/Zc^T, the expectation values in Eq.(p5D are short hand 
notations which express the following traces 

{H) = Tr[DH] , S = {-kBln{D)) = Tr[-kBDln{D)] , {N) = Tr[DN] . (36) 

Here ks is the Boltzmann constant, T the temperature, fi the baryon chemical potential, H 
the model Hamiltonian of the system, and N the baryonic particles number operator, given 
as usual by: 

N = [dx : ip-f^ifj : . (37) 



For a theory of free fields, i.e. when H reduces to a diagonal bilinear form in the field 
operators, computing the thermodynamics of the present system is a common textbook 
exercise |]T3. The partition function Z takes then a simple form which allows to express, for 
instance, the entropy S as Jl^ 

S= - kB^Y. [Up)lnUp} - (1 + /p(p))/n(l + f,{p))] 

p,p=iT,a 

-kB [fMinfM - (1 - fmn^ - fpim , m 

where /^^o- are Bose occupation numbers for the pion and the sigma modes, respectively, 
while fpp are Fermi occupation numbers for the baryons. For an interacting system of 
particles, the solution for the partition function Z as well as for the entropy S, as sketched 
above, persists in its simple and sympathetic form only in the case of an independent quasi- 
particle picture. The Hamiltonian, in this case, is split into a diagonal bilinear part which 
is fully admitted into the thermodynamics, and a residual interaction which is subsequently 
treated in a perturbation. There exist several ways of realizing a quasi-particle picture. 
These are based, for instance, on the self-consistent approximations a la Hartree-Bogoliubov 
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or a la Hartree-Fock-BogoliubovQ. Since we have established in the previous sections that, for 
symmetry requirements, the HFB basis is needed, we then chose to realize the independent 
quasi-particle picture in this latter. For this purpose, a thermal Bogoliubov rotation is 
applied to each field in the Hamiltonian according to: 

a+(T) = u^{p,T)a^ - v^{p,T)a.p ; /5+(r) = Ua{p,T)ht, - Vaip,T)b^p - wo{T)5{p) 



CpriT) = uf{p, T)cl^ - rvpip, T)d_ 



pr ) 



(T) = uf{p, T)d;^ + rvpip, T)c 



-pr 



(39) 



Here too the canonical normalization of the thermal amplitudes Un{p, T), Ua{p, T), Vt,{p, T), 
fo-(j9, T), uf{p,T) and vf{p,T) is assumed 

u^{p,Tf -v^{p,Tf = 1 , u^{p,Tf - v^{p,Tf = I , UFip,Tf + VFip,Tf = l. (40) 

Except for the temperature dependence, the above transformations strictly follow the zero 
temperature ones performed previously. Using Bloch - De Dominicis theorem |jl5|, it is 
straightforward to compute the traces (H) and (A^). The entropy 5* keeps its form in Eq.(^8D 
with the occupation factors given now in terms of the self-consistent thermal quasi-energies 
by: 



1 



A(g)= [exp {(3S^{q) 
/^(g)= [exp(/3(£:J(g)-/i)) + 1 
The thermal average of H reads: 



Uq)=[exp{(3Sj{q) 



1 



-1 



Mq)= exp(/3(fj(g)+/x)) + 1 



(41) 



+ 7 

3A2 
+ — 



(27r)3 2 
dq 
(2^ 



3 (1 + 2Uq)) K(g, TY + v^{q, Tf) + (1 + 2/<,(g)) (u^iq, Tf + v^. T) 



E, {uF{q,TY - VFiq,Ty] - 2g{a)TUFiq,T)vFiq,T) 



2/;i-\4 
T 



(42) 



With the definitions 



(ix (7ri(x)7rj(x)) 



(ix (o"(x)cr(x)) 



dx (^(x)^(x)) . (43) 



Minimizing the grand potential Q with respect to the shift parameter wo{T) gives an equation 
for the chiral condensate 



Sn ^ u^{0,T)+v,iO,T) 



S\'{a)T{lI + lI)-gI^ + \'m + f^'{a)r-c =0 (44) 



On the other hand, the variations of fl with respect to the thermal Bogoliubov amplitudes 
UTr{q,T), u„{q,T) and UF{q,T), while keeping the unitarity constraints Eq.(|40|) satisfied, 
lead to 



■^Indeed both of them deliver self-consistent operator bases which can keep the bilinear part of 
the Hamiltonian diagonal, thus allowing an independent quasi-particle picture. 
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UgU^{q,T)v^{q,T) + 



cliq) 



LJgUa{q,T)va{q,T) + 



A" {u^{q,T)+v^{T,q)y 
2 2uJn 



3A^ Mq,T)+v,{q,T)f 
2 2ujn 



5/J + /J + {arr 



H +1 



T 



{a)l 



5uF{q,T) 



cl{q) 



2EpUF{q, T) VF{q, T) + g{a)T Upiq, Tf - Vpiq, Tf 



(45) 



It is interesting to notice that the stationarity conditions hsted above ensure at the same 
time the independent particle picture for the finite temperature HFB solution, giving by 
this a consistency to the whole approach 



{[ai{q), [H,aj{p)]]) oc cl{q)6q,p6ij = 
mq),[H,P{p)]])cxcl{q)Sq,,^0 
{{Dr{q),{H, Cs{p)}}) oc 4(g)5,,/,, = 



(46) 



The thermal quasi-particle energies which correspond to the harmonic excitations of the 
system can be computed from 



{[a.iq),[H,anpm=£^{q)dqJ^^ 



u^{q,Tf + v^{q,Tf + 



mq)AH:P^{pm = £Uq)S, 



u,{q,Tf + v,{q,Tf + 3A 



[un{q,T) +v^{q,T)f 
2uj„ 



2{u a{q,T) +v„{q,T)f 

2Ua 



{{Cr{q).{H,Ct{p)}}) = £l{q)5q,,5rs 

= Eq (upiq, Tf - vpiq, Tf) - 2g{a)TUF{q, T)vF{q, T) 



57j + 7j + {a)l 
(47) 



Comparing the above results with those obtained earlier for the zero temperature case, one 
can see that the equations describing both stationarity conditions, namely the ground state 
energy for the first and the thermal equilibrium for the second, have kept the same form, 
apart of course from the presence of the thermal occupation factors in the latter. One can 
also check that at the minimum of Q, the quasi-particle energies read 



^J,<7(?) =^gK,a(g,T) -i;^,^(g,T)]^ , Sl{q) = ^ q^ + g'^ {0^)1 ^ 
which lead as well to the familiar expressions of the thermal tadpoles 



A 1 + 2/^, .(g) 
(27r)3 2 SIM ' 



ll = 27^^(0) 



d\ friq) + hiq) - 1 



{2ny 



2SM 



(48) 



(49) 



14 



Finally the finite temperature BCS solutions take the form of the following four coupled and 
self-consistent equations: 











{a)T 


c 




= -9{^)t ■ 



3A2(a) 



T tT 



li + 



3A2(a) 



T tT 



li + 



T / 



(50) 



As in the zero temperature case, the finite temperature BCS equations lead to a dy- 
namical mass generation and thus to a violation of the Goldstone theorem. Therefore, one 
needs to go one step further and consider those necessary thermal fluctuations present in the 
residual interaction. One way of bringing about these effects consists in taking the limit of 
a weakly interacting system which leads to a linearization of the thermal TDHF equations. 



Such a solution is known to correspond to the thermal RPA approximation (see ref. |16 
for details). The TRPA equation are very similar in form to the zero temperature Rowe 
equations of motion Eq.(^). The expectation values are, however, not realized anymore on 
the RPA ground state. They correspond, instead, to traces taken on the grand canonical 
ensemble, such that 



( H,Qt ) 



(51) 



The evaluation of such traces is a task which again can be handled very well by means of 
the Bloch-De Dominicis theorem. The net result is the following eigenvalue problem 



A^{quq2) S^(gi,g2) 



uj,{T)N\q- 



V.(gi,T) 



(52) 

where and are this time 4x4 matrices which have only an implicit dependence on 
temperature inherited from the temperature dependence of the thermal HFB basis. They 
are given by 



/^J(o) 









£^{qi)+£^Aqi) 


















S{qi - 52) + B^{qi,q2) 



(53) 



The effect of the residual interaction which is responsible for the RPA re-scattering is encoded 
in the matrix which takes the form 



Q2) 



2£J(0) 



/ 2A2(a)T7^^(g2) 

2X^a)Tn^{qi) 2X'n^{q,)n^{q2) 

2A2(a)T7^(gl) 2X' {q,)n^ (q^) 

V gi2n)-l 



2X^{a)Tn^iq2) 
2X'n^iq,)n^iq2] 
2X^n^{q,)n^{q2] 




7^(27r) 






-f \ 



(54) 
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Here TZ^iq) has an implicit temperature dependence and is of the same form as in Eq.(^). 

The exphcit temperature dependence of the RPA equations is carried by the norm matrix 
N"^ . This as well as the RPA amplitudes Uy{q,T) and Vy{q,T) are given by: 



Nj{q) 
-Nj{q) 



XPiq,T) 
Xi'Kq^T) 
\Xjf\q,T)J 



VMT) 



/Yy\q,T)\ 

Y('){q,T) 
\Y(^Kq,T)J 

(55) 



Now NJ is a 4 X 4 diagonal matrix, having as diagonal elements: 

Nliq)u = 1 , 

Njiq)22 = l + fAQ) + Uq) , 

Nl{q)33 = fAq)-Uq) , 

Njiq),, = 1 - Mq) - Mq) . 



(56) 



The difference between Nj and Id, the analog matrix at zero temperature Eq. (|27D , repre- 
sents the single formal departure of the thermal RPA equations from the zero temperature 
ones. Finally one can proceed to the resolution of the RPA eigenvalue problem. The RPA 
frequencies, as in the zero temperature case, can be read from the characteristic equation 



utiT) = Si ^(0) + g'^',,{ut{T)) + 



1 - 2X^J:Uu^I{T)) 



(57) 



where 



1 



1 - fFip) - fpip) 



dp 



1 + Up) + Up) 



£l{p)-£l{p) Up) -Up) 



+ 



'2snp)^^ip) ^1{T)-{£^{p)-£^Ap))' 

Now, like in the zero temperature case and making use of the identities 



KM 



one gets, after inspection of the gap equations, the following 



W^) 



(58) 



(59) 



^1{T) 



+ 9' 



T 



'FF 



S^f(O) 



2A2 



+ 



£P{0)-£r{0)] fe(^^(T))-SL(0) 



2X'j:U^iiT)) 



(60) 



Equation ( pOD exhibits clearly a zero energy solution for all temperatures below the transition 
to the Wigner Weyl phase. At the transition and beyond ((a) = 0), this solution is not 
normalizable, as can be easily seen from the norm matrix N"^. 



16 



Before closing this paper, it is certainly worth briefly showing how robust is the field- 
to-current mapping in systematically preserving the symmetry. A situation of particular 
interest in the nuclear problem is the renormalization of the fermionic four-point function 
by means of the so-called short range correlation. This can be addressed via a chiral invariant 
fermionic quartic interaction of Nambu-Jona-Lasinio (NJL) type 



+ ('0«75'r'0) 



(61) 



where the NJL-coupling G plays the role of the g' Migdal-parameter. The BCS gap equations 
are now further coupled by the insertion of a fermionic tadpole contribution to the fermion 
quasi-mass such that 





= /x^ + 




= /x2 + 3A2 


{a)T 


c 

-^^ 




= -9{^)t - 



5/J + /J 



2^F 



(62) 



where G — — 2G(1 + ^). Here again, the asymptotic Goldstone pion is generated by the 
previously introduced parity-mixing RPA. This time one gets the following 



+ 



9' 






1 - GEl^iuliT)) 


1 - GE^^(O); 



2A2 



+ 



gT2i 



(o)-^^r(o)i fek^m)-sL(o) 



1 - 2X^T.Uul{T)) 



(63) 



Figures (3) give the diagrammatic representation of the HFB-RPA dynamics in this case. Of 
course, further extensions of the model to other chirally invariant boson-fermion couplings 
are possible. 



IV. CONCLUSION 

In conclusion, we have presented an extension of the ficld-to-current mapping, introduced 
earlier, to a baryon rich regime at finite temperature. The mapping was made dynamical 
in the Hartree-Fock-Bogohubov-RPA approximation. The Goldstone theorem was fulfilled 
exactly, although the dynamics was not sorted-out according to neither an expansion in the 
available coupling constants (A, G) nor in the arbitrary charge and fiavor Np numbers. 
In fact, it can be easily appreciated, from the four coupled BCS equations for instance, that 
this approach is highly order mixing. 

All over this work, we did not consider the p-wave renormalization of the pion since it is 
well known that it does not endanger the Goldstone nature of the pion. Instead, we fixed 
our attention on the subtle situation of the s-wave renormalization where potential artefacts, 
linked to the dynamical mass generation, may alter the Goldstone character of the pion. This 
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was indeed visible from the BCS solution in which all the quasi-particle states were massive. 
Therefore it is certainly legitimate to state here that the HFB ground-state is not a viable 
vacuum and should not be tolerated as an approximate ground-state for nuclear matter 
saturation studies for instance. This point, although conceptually crucial and certainly 
quantitatively important, is usually not observed in the literature. As an alternative, we 
have presented the RPA ground-state, as a qualitatively correct approximate ground-state, 
in which the Nambu-Goldstone phase is exactly realized. On the other hand, we have also 
indicated that this approach can support whatever chirally invariant refinement made to 
the interaction Lagrangian. We hope that this will give a solid platform for a program of 
a quantitative assessment of the nuclear equation of state, beyond the usually considered 
classical mean-field calculations in QHD models. It should be also mentioned that the 
present work is of some relevance to the undergoing quantum mean-field calculations, based 
on chirally-invariant contact nuclear-forces, undertaken by Nikolaus, Hoch and Madland 
P!7| . According to these authors, the inclusion of the pion dynamical contributions is a step 
which one needs to envisage seriously in the future [|1^. The non-linear chiral realization a 



la heavy baryon ChPT is certainly not well suited for a non-perturbative treatment of the 
dynamics. Therefore we see in the present prescription an interesting framework for such a 
program. 
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FIGURES 




= .1 

FIG. 1. Diagrammatic representation of the four coupled BCS solutions. The dashed line 
denotes the self-consistent quasi-pion propagator, the solid line the quasi-sigma, the double solid 
line the quasi-fermion, and the wavy line the two point Green's function of the bare field a of the 
Lagrangian density in Eq.(||). 
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FIG. 2. Upper part: The Dyson equation for the physical pion (thick dashed Unes) for which 
the mass operator has been extracted from the scattering of the quasi-particles in a RPA equation. 
Lower part: The scattering equation for a pair of quasi-sigma (thin full lines) and quasi-pion (thin 
dashed lines). 
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FIG. 3. The full class of diagrams for the asymptotic Goldstone pion in the context of the 
extension discussed in the text. 
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